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Our work in directional solidification has been in the following arenas: 

! Dynamics of dendrites including rigorous mathematical analysis of the 
resulting equations. 

2. Examination of the near-structurally unstable features of the mathemati- 
cally related Hele-Shaw dynamics 

3. Numerical studies of steady temperature distribution in a vertical Bndg- 
man device. 


4 . 

5. 


Numerical study of transient effects in a 


vertical Bridgman device. 


Asymptotic treatment of quasi-steady operation of a vertical Bridgman 
furnace for large Rayleigh numbers and small Biot number m 3 


6. Understanding of Mullins-Sererka transition in a Bridgman device when 
fluid dynamics is accounted for. 


1 ' overlhe pS^upport, we have developed and refined a new technique to 
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in the complex plane is almost completely undeveloped- not much is known ex- 
cept the classic results of Cauchy-Kowaleski for first order equations. However, 
in the problems we encountered, the formal inner-outer asymptotic matching 
procedure sometimes requires us to have theoretical knowledge about existence 
of solutions of such parameter-free higher order complex partial differential equa- 
tion. Formal scaling arguments that give rise to important physical results on 
interface dynamics, assume the existence of such solutions. We have now made 
serious progress to extend the general theory of complex partial differential 
equation in a fundamental way [4] that justifies and extends the formal results 
obtained previously in [l]-[3]. Earlier, we [5] had rigorously obtained results 
on singularities of similarity solutions. Together with existence and uniqueness 
results in [4], we now have a generic information on the location of complex 
singularities for a wide-variety of situations. 

2. Examination of the near-structurally unstable aspects 

Usually, in a physical problem, small terms of the equation are assumed to 
have small effects on the solution. When this fails, the problem is deemed to be 
structurally unstable. A structurally unstable system is of course unphysical, 
since a physically realistic system always involve small effects that are neglected 
in the model. What is not so clear a priori is what might be expected for a 
system that is not structurally unstable, but close to it. It has been argued in 
a recent review by one of us [6] that both viscous fingering for small surface 
tension, as well as! dendrites for small anisotropy and small surface energy are 
examples of near structurally unstable system. For such systems, we have shown 
small terms in the equation that would seem reasonable to ignore can sometimes 
play a disproportionately large influence on the solution. Many surprising re- 
sults in viscous fingering and dendrites have been interpreted in this way. In 
particular it is argued that convection and many other effects that are estimated 
to be small in any carefully designed dendrite experiment can be unexpectedly 
important, despite the estimates. 

3. Numerical studies of steady Bridgman problem 

We have continued further extending work [7] in this direction. Attention 
has been placed on resolving the influence of thin ampoule walls on the steady- 
growth of crystals in a vertical Bridgman device. An asymptotic theory has 
been developed and tested against direct numerical simulation. There are three 
cases to consider. 

We assume that in the regions of the ovens, a Newton law of heat transfer 
applies at the outer wall of the ampoule. We then verify that this law can be 
applied to the inner wall when the heat conductivities of the ampoule and the 
melt or crystal are of comparable magnitude. In other words the ampoule is 

transparent to the heat flow conditions. 

For the region of the ampoule near the interface between a heater zone 
and the insulated region, there is an abrupt chnage in the heat transfer from 
a Newton law of heat transfer to no conduction. The thin ampoule has the 
effect of smearing the heat transfer coefficient so that a modified Newton law of 
cooling applies at the inner wall. An approximate version of this law was used 
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in our previous studies of steady crystal growth. The implication of the results 
is that thick ampoule walls lead to a degradation in radial segregation. 

The third and final case has proved more problematic. This case occurs 
where the crystal interface meets the ampoule wall. Previous numerical studies 
have indictated that heat conduction in the ampoule causes the interface to 
curve near the wall. What is overlooked is the singular nature of the interface at 
the wall We have found a similarity solution, but it does not meet the matching 
requirements for the temperature away from interface-wall region. Currently, 
a numerical study is being employed to identify the nature of the singu arity. 
While finite element methods are often capable of treating singularities in the 
solution domain, it is possible that there are serious errors in previous numerical 
treatments by not taking the singular nature of the interface’s contact with the 
ampoule in an appropriate way. 

Our plan is to resolve the singular nature of the interface before completing 
three-dimensional studies. 


4. Unsteady development of vertical Bridgman device 

A highly efficient code has been used to study the transient behavior of the 
temperature distribution in the Bridgman device [5]. In particular, we have 
identified the requirements on the length of the ampoule in order to guarantee a 
phase of steady crystal growth. Direct comparison with popular one-dimensional 
models reveal their short-comings. 

Some improvements have been made to the code: an exphct Adams-Bashfor 1 
method is now applied to both the convective terms and the interface motion, 
while Crank-Nicolson is still used to treat the rapid diffusion of heat. Currently 
tests are being conducted to to see whether direct solvers can replace the ADI 
splitting currently used in the code. If successful, and the prospects look good 
at the moment, direct solvers can then be applied for full three-dimensional 

simulations. . _ „ . , , . 

5. Asymptotic treatment of 3-D Bridgman with fluid motion 

We have extended the previous results reported by Tanveer [’94, Phys. Fluids 
A1 and Foster [8] for axi-symmetric situations into two ways. First, we have 
made the analysis fully three-dimensional, allowing for ampoule about which 
there is an azimuthal variation in the heat transfer through the wide walls, or 
azimuthal variation due to a slight tilt of the ampoule axis [9] We find that 
the three-dimensional effects in every case increase the horizontal segragation 
in the crystal. We also determined that the optimization criteria found first 
bv Tanveer extend in verv similar form to the thee-dimensional situation, but 
now involve zeros of the >! Bessel function rather than J 0 . Secondly, we haw 
explored effects of rapid rotation of the ampoule about the vertical axis [10]. 
We find, in agreement with other recent numerical work, that in fact rotation 
of the ampoule worsens the radial segregation in the crystal. W hat effect that 
rotation has on the onset of morphological instability awaits further work. 


6. Mullins-Sererka transition in Bridgman device 
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St"at obtained in an unmstr icted — ^ £%££ 
obtained by a spectral scheme. 
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